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Abstract 

A stochastic representation for the solutions of the Poisson-Vlasov 
equation is obtained. The representation involves both an exponential 
and a branching process. The stochastic representation, besides pro- 
viding an alternative existence proof and an intuitive characterization 
of the solutions, may also be used to obtain an intrinsic definition of 
the fluctuations. 



1 Introduction 

The solutions of linear elliptic and parabolic equations, both with Cauchy 
and Dirichlet boundary conditions, have a probabilistic interpretation, which 
not only provides intuition on the nature of the problems described by these 
equations, but is also quite useful in the proof of general theorems. This 
is a very classical field which may be traced back to the work of Courant, 
Friedrichs and Lewy [1] in the 20 's. In spite of the pioneering work of McKean 
[2], the question of whether useful probabilistic representations could also be 
found for a large class of nonlinear equations remained an essentially open 
problem for many years. 

It was only in the 90's that, with the work of Dynkin[3] such a theory 
started to take shape. For nonlinear diffusion processes, the branching exit 
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Markov systems, that is, processes that involve diffusion and branching, seem 
to play the same role as Brownian motion in the linear equations. However 
the theory is still limited to some classes of nonlinearities and there is much 
room for further mathematical improvement. 

Another field, where considerable recent advances were achieved, was the 
probabilistic representation of the Fourier transformed Navier-Stokes equa- 
tion, first with the work of LeJan and Sznitman[5], later followed by extensive 
developments of the Oregon school [6J ^ In all cases the stochastic rep- 
resentation defines a process for which the mean values of some functionals 
coincide with the solution of the deterministic equation. 

Stochastic representations, in addition to its intrinsic mathematical rele- 
vance, have several practical implications: 

(i) They provide an intuitive characterization of the equation solutions; 

(ii) They provide a calculation tool which may replace, for example, the 
need for very fine integration grids at high Reynolds numbers; 

(iii) By associating a stochastic process to the solutions of the equation, 
they provide an intrinsic characterization of the nature of the fluctuations 
associated to the physical system. In some cases the stochastic process is 
essentially unique, in others there is a class of processes with means leading to 
the same solution. The physical significance of this feature is worth exploring. 

A field where stochastic representations have not yet been developed (and 
where for the practical applications cited above they might be useful) is the 
field of kinetic equations for charged fluids. As a first step towards this 
goal, a stochastic representation is here constructed for the solutions of the 
Poisson-Vlasov equation. 

The comments in the final section point towards future work, in partic- 
ular on how a stochastic representation may be used for a characterization 
of fluctuations, alternative to existing methods. This is what we call the 
stochastic principle. 

2 Stochastic representation and existence 

Consider a Poisson-Vlasov equation in 3+1 space-time dimensions 

df ^ e 

^+ V -v./ - -V.$ ■ V./ = (1) 
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with 



A,*/ = -47r |e j f (x, V, £v - eps (x) | (2) 

Pb (^x^ being a background charge density. 
Passing to the Fourier transform 

F{U) = -^ [ d'vf{v,t)S-^ (3) 

[271) J 

with 7]= (^x^v^ and ^ = ^^1,^2^ = (?i5'^2), one obtains 

= ^^^-S-V,,F(e,t) (4) 
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(^1) being the Fourier transform of pb (x). Changing variables to 

r = limt (5) 
where 7(|.^2|) is a positive continuous function satisfying 

7(161) = 1 if 161 <1 

7(161) > 161 if 161 >i 

leads to 

9F{^,t) 6 ^ , 47re 



dr 7(161) 



V,,F (6 d%F (6 - 6, r) 



with 6= ]|7y- Eq.Q written in integral form, is 

F(6r) = e"^-^«2F(6,6,0)-^ rc/se^""'^^-^«2 (7) 
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For convenience, a stochastic representation is going to be written for the 
following function 



(8) 



with A a constant and h{^i) a positive function to be specified later on. The 
integral equation for x (^ii ^2, t") is 
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6-^1 

7(161) 



7(161)' 

1 Pb 
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with 



and 



{\iirh*h) 



-1 



(10) 
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{\Ci\-'h*h) 

Eq.Q has a stochastic interpretation as an exponential process (with a 
time shift in the second variable) plus a branching process. P d^^i 
is the probability that, given a ^1 mode, one obtains a (^1 — branch- 
ing with ^[ in the volume (^1,^1 + d^^'i)- x{CijC2,t) is computed from the 
expectation value of a multiplicative functional associated to the processes. 
Convergence of the multiplicative functional hinges on the fulfilling of the 
following conditions 
>(6.g2.o)' 



(A) < 1 

(c) {\^i\-'h*h) </i(ei) 

Condition (C) is satisfied, for example, for 



(9) 
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Indeed computing ^ h* K) one obtains 



/•oo 

— {\ii\-^h*h)=ATic dr 
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(13) 

This integral is bounded by a constant for all therefore, choosing c 
sufficiently small, condition (C) is satisfied. 

Once h (^i) consistent with (C) is found, conditions (A) and (B) only put 
restrictions on the initial conditions and the background charge. Now one 
constructs the stochastic process X (i^i, ^2, t)- 

Because e"^"^ is the survival probability during time r of an exponential 
process with parameter A and \e~'^^ds the decay probability in the interval 
(s,s + (is), x(G)'?2i''") in Eq.Q is obtained as the expectation value of a 
multiplicative functional for the following backward-in-time process: 

Starting at (^i,^,^), a particle lives for an exponentially distributed 
time s up to time r — s. At its death a coin Is (probabilities \,\) is tossed. 
If /s = two new particles are born at time t — s with Fourier modes 

(ii - Cl, 6 + and (^^ , O) with probability density p (^1, Q .111^ = 1 

only the — i^i, (^2 + -^^5^^^) particle is born and the process also samples 

the background charge at Pb {^ij- Each one of the newborn particles contin- 
ues its backward-in-time evolution, following the same death and birth laws. 
When one of the particles of this tree reaches time zero it samples the initial 
condition. The multiplicative functional of the process is the product of the 
following contributions: 

- At each branching point where two particles are born, the coupling 
constant is 

- When only one particle is born and the process samples the background 
charge, the coupling is 



mX 
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- When one particle reaches time zero and samples the initial condition 
the coupling is 

50(6,6) = ^^^ (16) 

The multiplicative functional is the product of all these couplings for each 
realization of the process X (^1, ^2, ^), this process being obtained as the limit 
of the following iterative process 

Then, x (Ci) ^2) ''') is the expectation value of the functional. 

X(6,6,r) =E|n(5o^;---) (^i^;---) {g29'2---)] (17) 

For example, for the realization in Fig.l the contribution to the multi- 
plicative functional is 

(6 - - 'Ci , ^2) 9q (Cl - Ci') ^3) 9o (Ci , 0) 



and 



ki — k-\- si 



7(161) ^ 
A;2 = A;i + (S2 - Si) 
^3 = (S3 - Si) 

With the conditions (A) and (B), choosing 

_ Sttc^ 
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and 

c < 
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Figure 1: A sample path of the stochastic process 



the absolute value of all coupling constants is bounded by one. The branch- 
ing process, being identical to a Galton- Watson process, terminates with 
probability one and the number of inputs to the functional is finite (with 
probability one). With the bounds on the coupling constants, the multi- 
plicative functional is bounded by one in absolute value almost surely. 

Once a stochastic representation is obtained for x {^i, ^2, t), one also has, 
by dS]), a stochastic representation for the solution of the Fourier-transformed 
Poisson-Vlasov equation. The results are summarized in the following : 

Theorem 2.1 - There is a stochastic representation for the Fourier- 
transformed solution of the Poisson-Vlasov equation F{C,i,^2,t) for any ar- 
bitrary finite value of the arguments, provided the initial conditions at time 
zero and the background charge satisfy the boundedness conditions (A) and 
(B). 
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As a corollary one also infers an existence result for (arbitrarily large) 
finite time. Notice that existence by the stochastic representation method re- 
quires only boundedness conditions on the initial conditions and background 
charge and not any strict smoothness properties. 

3 Fluctuations and the stochastic principle. 
A comment 

In the past, the fluctuation spectrum of charged fluids was studied either by 
the BBGKY hierarchy derived from the Liouville or Klimontovich equations, 
with some sort of closure approximation, or by direct approximations to the 
N-body partition function or by models of dressed test particles, etc. (see 
reviews in |j9j [IQj). Alternatively, by linearizing the Vlasov equation about 
a stable solution and diagonalizing the Hamiltonian, a canonical partition 
function may be used to compute correlation functions [TT] . 

However, one should remember that, as a model for charged fluids, the 
Vlasov equation is just a mean-field collisionless theory. Therefore, it is 
unlikely that, by itself, it will contain full information on the fluctuation 
spectrum. Kinetic and fluid equations are obtained from the full particle 
dynamics in the 6N-dimensional phase-space by a chain of reductions. Along 
the way, information on the actual nature of fluctuations and turbulence 
may have been lost. An accurate model of turbulence may exist at some 
intermediate (mesoscopic) level, but not necessarily in the final mean-field 
equation. 

When a stochastic representation is constructed, one obtains a process 
for which the mean value is the solution of the mean-field equation. The 
process itself contains more information. This does not mean, of course, that 
the process is an accurate mesoscopic model of Nature, because we might be 
climbing up a path different from the one that led us down from the particle 
dynamics. 

Nevertheless, insofar as the stochastic representation is qualitatively unique 
and related to some reasonable iterative proces^, it provides a surrogate 
mesoscopic model from which fluctuations are easily computed. This is what 
we refer to as the stochastic principle. At the minimum, one might say that 

^Representations as those constructed for the Navier-Stokes equation and the one in 
this paper may be looked at as a stochastic version of Picard iteration 
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the stochastic principle provides another closure procedure. 
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